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Abstract 

A model of spin torque induced magnetization dynamics based upon semi-classical spin diffusion 
theory for a single layer nanocontact is presented. The model incorporates effects due to the 
current induced Oersted field and predicts the generation of a variety of spatially dependent, 
coherent, precessional magnetic wave structures. Directionally controllable collimated spin wave 
beams, vortex spiral waves, and localized standing waves are found to be excited by the interplay 
of the Oersted field and the orientation of an applied field. These fields act as a spin wave "corral" 
around the nanocontact that controls the propagation of spin waves in certain directions. 

PACS numbers: 75.40. Gb 85.75.-d 75.40.Mg 76.50.+g 75.30.Ds 75.75.+a 75.70.Cn 72.25.Ba 
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The flow of sufficiently large dc current through a thin, nanomagnetic structure can give 
rise to precessional magnetization dynamics at gigahertz frequencies This remarkable 
effect has attracted broad interest, both from the standpoint of fundamental physics and in 
the context of applications. The underlying physics of these spin torque devices is based 
upon the ability of thin ferromagnetic layers to act as spin filters when current flows through 
the layers. For spin torque effects to manifest, a source of spin polarized carriers with a 
component perpendicular to the magnetization of a layer is required. A typical spin torque 
multilayer has two primary magnetic layers: a fixed layer that acts as a spin "sieve" that 
induces a spin accumulation in a non-magnetic spacer layer, and an active layer that can 
respond dynamically when it absorbs the angular momentum from the accumulated spins. 
If the active layer has a uniform magnetization, a torque is produced only when the two 
layers are misaligned 

On the other hand, if the magnetization is not uniform, theory for even a single magnetic 
layer predicts a non-zero torque with resulting dynamics ^. Single layer spin torque the- 
ory was used to explain differential resistance data in mechanical nanocontact experiments 
5] and in lithographically defined nanopillars Q]. Theoretical studies have considered single 
layer nanocontact devices [3], but have not addressed the response of a physically realistic, 
finite sized nanocontact with its accompanying Oersted field generated by the dc current 
flowing through the device. 

In general, the total spin accumulation, and hence spin torque, in a magnetic thin film 
device arises from the lateral and longitudinal diffusion of spins, transverse and parallel to 
the current direction, respectively 4]- The oft used Slonczewski model of spin torque in 
trilayers, however, assumes a uniform spin accumulation and incorporates only longitudinal 
spin diffusion effects 

In this paper, we report a realistic treatment of a single layer device including lateral 
and longitudinal spin diffusion, the Oersted field, and a large enough sample geometry to 
capture novel behavior. Using a novel micromagnetic simulator, we demonstrate unexpected 
features of the response, including localized standing waves, vortex spiral waves, and, most 
strikingly, a weakly diffracting collimated beam of spin waves, the direction of which can be 
steered by changing the direction of an applied magnetic field. The formation of the beam 
appears to be a novel physical mechanism involving the hybridization of a localized standing 
wave and a vortex spiral wave. It has been previously shown that spin waves emitted from 
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FIG. 1: Single layer nanocontact device schematic. The spatially nonuniform Oersted field hoe is 
generated by the dc current /. "Py" = Permalloy (Ni8oFe2o)- 



nanocontact devices can be used to phase-lock two spin torque oscillators [8| . The ability to 
steer a spin wave beam with magnetic field could offer a method to selectively control phase 
locking of multiple spin torque oscillators in an array structure. 

The outline of this paper is as follows. First, we describe a two-dimensional (2D) model of 
spin torque in single layer, nanocontact devices. We then present micromagnetic simulations 
that demonstrate the wide variety of responses mentioned above. Finally, we explain the 
results of the simulations using a local formulation of the linear spin wave dispersion relation 
above and below the nanocontact, showing how the applied and Oersted fields act as a spin 
wave "corral" enabling directional control of anisotropic spin wave propagation and localized 
excitations. We also show that this corralling effect is limited neither to this particular model 
nor the details of spin torque. 

The physical system we analyze is pictured in Fig. [Hand is similar to the one theoretically 
studied in {4] except that we explicitly treat the finite contact area and Oersted fields. A 
single ferromagnetic Ni8oFe2o (Permalloy, Py) layer is adjacent to two copper (Cu) leads and 
an insulator. The current flows uniformly in the —z direction (electron flow is in the +z 
direction) from a right reservoir located at z = + through the cylindrical Cu lead of 
radius r^., across the Py layer of thickness t at 2; = 0, over to a left reservoir at a distance 
away from the magnetic layer [z = —ti)- The length t/j is an effective distance over which 
the current is assumed to maintain quasi-unidirectional flow. The magnetic layer is assumed 
to have infinite extent in the xy directions. When we refer to the region above (below) the 
nanocontact, we mean the positive (negative) direction along the y-axis. 

We calculate the spin accumulation due to current flow through a ferromagnet using the 
same method described in Ref. jj]. We consider the behavior of the spin accumulation in 
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the non-magnet in response to a non-uniform magnetization M = M^u, where Ms is the 



saturation magnetization. As in Ref. {4], we treat the case of small deviations u± away 
from the spatially averaged "equilibrium" direction u\\. The transverse and longitudinal 
components of the spin accumulation m = mj_ + my, where my points in the longitudinal 
direction of the steady state spin accumulation in the absence of any magnetic inhomogeneity, 
can be decoupled in this limit. We solve the multipoint boundary value longitudinal problem 
to find the longitudinal spin accumulation my and spin current Q^z for each interface. For 
deviations at a particular transverse wave vector [k^, ky), the transverse spin accumulation 
is then given in terms of the deviations in the magnetization and the longitudinal solution 

=^{u±}- " , 

Dk coth[l' K) + Wq ^-j^-^ 

^ = ikl + kl + l/ll,Y/^ 

where J^im^} and J^{u±} are the 2D Fourier transforms of the transverse spin accumulation 
and the magnetization transverse to the average, respectively (+ for right z = t interface 
and — for left z = interface). The decoupling of the longitudinal spin accumulation 
from the transverse spin accumulation is strictly valid only in the limit of small deviations 
from a uniform magnetization distribution. In our case, the deviations from uniformity are 
not small, so this treatment should be considered to be a first order approximation. The 
distance to the reservoir /' is or t/j, for the left and right interfaces respectively. The 
spin diffusion length Z^/ and diffusion constant D are material parameters for Cu; Wq is the 
effective interfacial spin absorption rate. 

To calculate the transverse spin accumulation in real, polar coordinate (r, 0) space, we 
take the inverse Fourier transform of eq. ([T]). Then the expression for m±, eq. ( lA.ll) (see 
appendix), is the convolution of the magnetization with a weakly singular kernel over the 
point contact region. This nonlocal formulation of the transverse spin accumulation can be 
interpreted as the lateral diffusion of spins interacting with a given magnetization distribu- 
tion. 

We find the quasi-steady-state spin accumulation for a given instantaneous magnetization 
distribution in the Py layer. This is justified because the ratio of the time scales for the 
diffusion of electrons to a steady state and for the magnetization dynamics is about 0.001. 
By formulating the calculation of the inhomogeneous transverse spin accumulation in terms 
of a simple convolution operation, we have greatly improved the speed of simulating this 



effect as compared to directly calculating the coupled magnetization and spin accumulation 

We use the average magnetization direction over the contact u\\, eq. (lA.2p . as the orien- 
tation of the longitudinal spin accumulation. Physically, this corresponds to the situation 
where a spin scattered a large number of times from the interface effectively "sees" the av- 
erage magnetization. The total spin accumulation inside the nanocontact m*, eq. (]A.3p . is 



then calculated by summing the contributions to the longitudinal and transverse accumu- 
lation from each interface. Because an insulator surrounds the Cu leads in our model, the 
spin accumulation outside the nanocontact is zero. 
The dynamical equation for the magnetization is 

dii -* -* 

— = —u X — OiU X (m X /icflf) + au X {u x m*), 



ics = ho - - g{r)(f) + r]V^u, a ° 



This is a modified 2D Landau-Lifshitz equation in dimensionless form with time normalized 
by 7/ioMs (7 is the gyromagnetic ratio, /io the permeability of free space). In eq. ([2]), space is 
normalized by r^,, fields and magnetization are normalized by Mg, a is the damping constant, 
rj = Dcx/i'yp-oMshrl) is the coefficient of the exchange term (-Dex is the exchange param- 
eter, h is Planck's constant divided by 2n), ho = /io[sin(6'o) cos(?/'o)5 sin(6'o) sin(?/'o)5 cos(6'o)] 
represents the canted, normalized applied field (^0 is measured from the positive ^-axis, ipo 
is measured counterclockwise from the positive x-axis), —UzZ is the axial demagnetizing 
term, and h j{r, 0) = —g{r)(j) is the nonuniform Oersted field due to the current density J(r) 
defined below. The driving torque is similar to the Slonczewski torque for a trilayer device 
[2! except that here, the torque is nonlocal. We have only considered the leading order axial 
component of the dipole field, —UzZ that is independent of film thickness t. Calculation of 
higher order terms involves integration over the ferromagnetic volume that greatly compli- 
cates the solution of eq. ([2]). Local approximations to higher order correction terms for such 
dipole fields have been derived by Arias and Mills [l^ . Such correction terms scale as t V ■ u 
and have quantitative significance for large film thicknesses. 

The Oersted field is significant in the nanocontact geometry, with a maximum magnitude 
on the order of 80 kA/m (1000 Oe) for r* = 40 nm and / = 21 niA. Thus, to ignore the 



Oersted field, as was done in previously presented multilayer simulations ll|, is not an 



appropriate approximation. We will demonstrate that the Oersted field significantly affects 



the response of the system. Now we derive the Oersted field for the geometry shown in Fi^ 

m 

The current density J{r', z') in dimensional units (r', z') is modeled as 
J(r', z') = 

f-^F{r') [H{z' - t - tR)H{t + tn + D,- z')- 



(3) 



H{-tL-z')H{z' + tL + D^)] 
- z^H{r')H{n - r')H{t + t« - z')H{z' + t^), 

J r'/r, < r' < J x < 

F(r ) = <j , H[x) = < 

I r*/r' r^, < r' I 1 a; > 

Two infinite parallel conductor plates of thickness Dp, separated a distance + t + t^, are 
connected by a wire of radius r^,. Current fiows into the wire from the right plate, through 
the magnetic layer via a cylinder modeling the point contact, and out of the wire into the 
left plate. The current fiow in the plates is assumed to be in the radial direction only. The 
coefficient of f in ([3]) models the magnitude of the current density in the plates, assumed 
uniform in the z direction. Outside of the region where the wire connects to the plates 
(r > r*), the current density falls off proportional to 1/r, whereas inside the region (r < r,,), 
the current is proportional to r. The current is assumed to flow uniformly down the wire 
modeled by the coefficient of i in ([3]). This model of current flow conserves the flux of 
current from the plates to the wire. We are interested in the value of the Oersted fleld Hj, 
due to the current density ([3]), in the center of the magnetic layer z = t/2. 

Using a vector potential representation along with Fourier and Hankel transforms, we 



solve for the Oersted fleld in dimensionless units (r, z) = (r'/r^,, z' /r^) 121]: 



infinite wire finite wire correction 



ly ff ry 

* [h{-) + h{-)]] (4) 



conductive plates 
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The integrals in (j3j) are 



Hp)- 



p/2 < p < 1 
l/2p 1 < p 



1 



h{p)= e-'^''/\osh{qz.)Mq)Mqp)-dq, 

Jo Q {•^j 

poo -I 

h{p)= / e-'^-'\l-e-^'')coMqz,)Uq)Uqp)-dq, 
Jo Q 

j-oa 1 

h{p)= / e-^'^/2(l-e-'^'^)cosh(gz,)J2(g)Ji(gp)-dg, 
Jo Q 

where a = (tji + t + tL)/r^, d = Dp/r^,, = {tji — t^)/2r^, are the normahzed wire length, 

conductive plate thickness, and the location of the magnetic layer, respectively. J„ are the 

nth-order Bessel functions of the first kind. Note that (jll), with I2 = I3 = I4 = 0, is the 

result for the magnetic field due to an infinitely long wire of radius r^, with no conductive 

plates. The magnitude of the Oersted field hj{r) = \hj{r, </))| in (jlj) depends only on r. 

By numerically evaluating the integrals ([5]), we show Mshj{r) for two contact sizes and 

a fixed current / = 30 mA in Fig. [2l The fields due to an infinitely long wire with no 

conductive plates are also shown in Fig. [2] as dashed curves. Equation (jlj) and the result 

for an infinite wire agree near the point contact but the decay of the fields is faster for the 

infinite wire case. 

/ = 30 mA, r» = 20 nm / = 30 mA, r, = 40 nm 
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FIG. 2: The magnitude of the Oersted magnetic field (jl]) (solid) and the Oersted field due to an 
infinitely long wire with no conductive plates (dashed). 



We implemented a numerical method to solve eq. (I2j) in polar coordinates, (see appendix). 
The calculations are rendered tractable by formulating the model in a nonuniform polar 
coordinate grid, allowing us to compute over a large domain (4.8 /im diameter disk) to avoid 
boundary spin wave reflections and with simulation times (3 ns) sufficiently long to ensure 



FIG. 3: Magnetization (uy) pattern and energy density of a thin film excited by dc current through 
a nanocontact with varied apphed field canting angle Oq = 1° , 10°, 21° in (a,d), (b,e), and (c,f), 
respectively. All other system parameters are fixed (/ = 29 mA, = 0). The circle in the 
center represents the boundary of the nanocontact. The domain is a square ten times the contact 
diameter per side. Magnitudes in each panel are normalized; positive values are yellow, negative 
values are black, with Uy oscillating between approximately ±0.8, ±0.6, ±0.65 in (a), (b), and (c), 
respectively. The peak energy density at r = 10 in (d), (e), and (f) is, in arbitrary units, 0.74, 1.0, 
and 0.15 respectively. Note the pinning of the vortex core in (a) and (d). The energy density plot 
in (e) clearly shows the weak diffraction of the spin wave beam. 

that we have determined the true steady state response. By evolving eq. in time with a 
nonuniform initial condition (where u is relaxed in the presence of the effective field only), we 
find that the magnetization settles into a quasiperiodic state due to the competition between 
the spin accumulation torque and the damping. All excitation frequencies are calculated 
from the time series of Uy averaged over the nanocontact by use of Fourier methods with a 
typical resolution of 0.75 GHz. We use the physical parameters listed in Table I of [4] except: 
ho = 1.1, = 40 nm, tR = 5 nm, ti = 75 nm, t = 15 nm, = 800 kA/m, D' = 0.001 mVs 
(diffusion rate in Py), If^^'^ = 5.5 nm (spin diffusion length in Py), and Dp = 50 nm. 

Our calculations show a variety of behaviors that depend on the physical parameters 
including vortex spiral waves (Figs. [3^, HJ:;), spin wave beams (Figs. [3]d, |4)d), and localized 
standing waves (Figs. [3t, S^). The direction of spin wave propagation in a canted field, 
^0 > 0, depends on the applied field azimuthal angle ipo where the spin wave angle ~ ipo — dO". 
The top panels depict the spatial variation of Uy at a specific time and the lower panels 
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FIG. 4: Magnetization (uy) pattern and energy density of a thin film excited by dc current through 
a nanocontact with varied dc current / = 18, 28.5, 36 mA in (a,d), (b,e), and (c,f), respectively. 
All other system parameters are fixed (^o = 10°, ■00 = 0°). Uy oscillates between approximately 
[-0.37,0.34], [-0.53,0.52], [-0.97,0.87] in (a), (b), and (c), respectively. The peak energy density 
at r = 10 in (d), (e), and (f) is, in arbitrary units, 0.0051, 0.35, and 1.0 respectively. 

show the energy density for each state. The energy density E[r, (p) is calculated by taking 
the time average of the squared magnitude of the magnetization transverse to the average 
direction E{r,(f)) = {\u x {u) p), where (/) = ^ JJ*^^ fdr, and r^, is the time by which the 
magnetization has settled into a steady, precessional state. 

Our choice of a relatively thick Py layer emphasizes the effect of the Oersted field. Larger 
currents are necessary to excite thicker layers, so that the Oersted fields are larger. Cal- 
culations with a thinner Py film still result in localized standing waves, spiral waves, and 
anisotropic waves, but the collimated beam is more difficult to excite. 

The particular response excited can be explained, in part, by appealing to the dispersion 
relation for a spin wave propagating in a canted uniform field, eq. (lA.4p . The local fields in 
the regions above and below the nanocontact (on either side of the nanocontact in the plane 
of the magnetic film) are of different magnitudes and orientations due to the presence of the 
nonuniform Oersted field —g{r)(j). The Oersted field acts as a "corral" and effectively lifts the 
spatial degeneracy of the dispersion relation immediately above and below the nanocontact 
so that, at a given frequency, the spin waves propagate in one direction and evanesce in the 
other. 

The response away from the nanocontact does not strongly depend on the details of 
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spin torque except that the spin torque is locahzed at the nanocontact. To show this, we 
solved eq. ([2]) with a locahzed ac applied field {h^dr, (j),t) = /lac sin(27r/acr)^, r < 0.15, 
elsewhere}, neglecting spin torque. The response and the associated dispersion curves are 
depicted in Figs. E^-c for / = 30 niA, 9q = 10° , ipQ = 0° , and /^ac = 1- We use the dispersion 
relation with the local fields evaluated at r = 1 to approximate which wavenumbers can 
propagate above (solid blue) and below (dashed black) the nanocontact. The dash-dotted 
red curve is the far field dispersion curve where the Oersted field is negligible. The Oersted 
field creates a gap between the dispersion curves above and below the nanocontact. The type 
of response excited depends on the driving frequency /ac and its relation to the ferromagnetic 
resonance (FMR) frequencies (eq. (lA.4p with no exchange contribution, r/ = 0) above (/j^mr)) 
below (/fmr); ^iid far away (/Imr) fro^i the nanocontact. When the driving frequency is 
between fpy^^i /fmr^ ^^e excitation is a standing wave (see Fig. [5^). When the driving 
frequency is between /|^r /fmr; excitation forms a spin wave beam (see Fig. [5)d). 
When the driving frequency is above /j)mr5 ^^e excitation is a nonlocalized propagating 
wave (see Fig. [5b). This is the "corral" effect; the applied field and the Oersted field act in 
concert to modify the availability of spin wave states in close proximity to the nanocontact, 
depending upon the excitation frequency. In the locations where the Oersted and applied 
fields add (above the nanocontact), the dispersion relation is shifted upward in frequency, 
thereby acting as a "fence" to block spin wave propagation if the excitation frequency is 
below /fmr- Where the Oersted and applied fields cancel, the dispersion relation is shifted 
downward in frequency, thereby acting as a trap, localizing the spin wave excitations until 
the excitation frequency exceeds /Imr' whereupon the trap is no longer operative, and the 
spin waves are guided downward. The mechanism leading to the standing wave is a linear 
phenomenon, in contradistinction to previous calculations where nonlinear effects lead to 
the formation of a localized magnetic excitation in a nanocontact geometry 13|. Even in 
the presence of a small driving field (e.g., /lac = 0.001), both the localized wave and the spin 
wave beam can be excited, precluding any strong role of nonlinearity in the localization of 
the response. 

We can interpret the different spatial responses in Figs. [3], H] by plotting the excitation 
spectrum as a function of applied field angle (Fig. [6]) and applied dc current (Fig. [7j). All 
future local FMR and dispersion curve calculations assume the Oersted field —g{r)(j) is 
evaluated at the radius r = tq = 1.5. For currents close to the excitation threshold, the 
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FIG. 5: Band structure, excited frequencies and response energy density for three different driving 
frequencies (a) /ac = 10 GHz, (b) /ac = 16 GHz, (c) /ac = 25 GHz. The fined/hohow circles 
correspond to the numerically determined wavenumber above/below the nanocontact and agree 
with the linear dispersion relation. The FMR frequencies above and below the nanocontact are 
/fmr ~ ^0.2 GHz and /j^mr = 5.6 GHz, respectively. When the drive frequency is in the gap 
between the upper and lower dispersion relations, the second harmonic 2/ac is also excited. 

response spectrum has a single frequency peak that corresponds to vortex spiral waves for 
small angles and a localized standing wave for larger angles (see the inset of Fig. [H]). The 
lowest available frequency response is excited. For small angles, the localized wave is not 
available because /Imr (red, dash-dotted curve) is less than f^uR (white, dashed curve) so 
spin waves can propagate to the far field and the vortex spiral wave is excited. Above the 
critical angle 9,^ = sin~^[(7(ro)/(2/io)] where /Imr = /fmr radius r = tq, the localized 

wave with a much lower frequency is excited because its frequency is below /fmr (see the 
vertical dotted lines in Fig. E]). 

For larger applied currents, nonlinear effects appear to be important. In Fig. El the vortex 
spiral wave is excited for small angles for the same reasons mentioned previously. However, 
the bifurcation in mode behavior occurs for angles larger than 6^, where, now, there are two 
distinct frequency peak branches in the spectrum that are not harmonically related. These 
branches lie just above /fmr /fmr (solid, blue curve). In the small amplitude case (Fig. 
EJinset), these branches correspond to the localized standing wave and the vortex spiral wave, 
respectively. Here, nonlinearity appears to spawn a hybridization of the standing wave and 
the vortex spiral wave resulting in the spin wave beam (see Figs. [3]d and|4)D). As the angle 
is increased, the power associated with each branch changes from being predominantly in 
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the vortex spiral wave branch to being predominantly in the localized wave branch and 
back again resulting in a visual change of the energy density as shown in Fig. ^ where the 
standing wave is dominating the response. 

This nonlinear hybridization is further investigated in Fig. [7] where the frequency spec- 
trum is plotted as a function of applied current (a) and associated local spin wave dispersion 
curves are plotted (b-d). For small currents, the localized standing wave is excited due to the 
corral effect (see Fig. [7)d and compare with Fig. [5^). For large currents, anisotropic waves 
are excited because the corral is partially open (see Fig. [711 ). For intermediate currents, 
between the vertical dashed lines, both the localized wave and the spiral wave are excited at 
the same time, but with different frequencies /" and It is only within this current range 
that we observe the spin wave beam. There is a single wavenumber associated with these 
two frequencies, k~ , if the dispersion curve below (above) the contact is associated with the 
localized standing wave (vortex spiral wave) (see Fig. [7t). The spin wave beam corresponds 
to a wavenumber associated with the dispersion curve below the contact but with the 
larger, vortex spiral wave frequency These two frequencies and two wavenumbers form 
a triad in the dispersion diagram [7b. We interpret this behavior in the following way: 



Given the coarse approximations made in this analysis (linear, ID behavior and pointwise 
evaluation of the Oersted field, etc.), the discussion provides a qualitative description of the 
response selection in single layer nanocontacts. 

The results predicted here for nanocontacts stand in stark contrast to calculations in- 
volving nanopillar structures [l^ where excitations are confined to a nanomagnet. In the 
latter case, the Oersted field was shown to effect a small perturbation to the magnetization 
dynamics. Here, the Oersted field gives rise to the "corral" effect and is fundamental to the 
formation of magnetic excitations. 

We have also performed numerical simulations in trilayer structures (ferromag- 
net/nonmagnet/ferromagnet) neglecting lateral diffusion and nonlocal effects and using the 
Slonczewski torque M in both our own simulator and in a conventional micromagnetic sim- 

n 

ulation package (OOMMF) by a colleague [15[. We find similar qualitative behavior such 
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FIG. 6: Nanocontact frequency spectrum as a function of applied field polar canting angle 
along with the local FMR frequencies above (/j?mr)' below (/pj^j^), and far away (/f^r) from the 
nanocontact. Main plot current: I = 29 mA; inset plot current: / = 18 mA. 

as the corralling effect and a spin wave beam. However, we do not observe stable vortex 
spiral waves, and the spatial dependence of the excitations are different. For spin torque 
magnitudes comparable to that of the single layer torque discussed in this work, an unphys- 
ical, large spin transfer efficiency of 0.7 was required. Theoretical calculations suggest that 
the spin transfer efficiency is between 0.25 and 0.4 . Thus, we find that the torque 

generated by a single layer with nonuniform magnetization can exceed that expected for a 
uniformly magnetized layer in a trilayer structure, suggesting that inclusion of spin diffusion 
effects considered here in trilayer calculations may have significant quantitative impact. We 
note that the semiclassical transport calculations in ^ used here do not require the use of 
the spin transfer efficiency, an ad hoc parameter, because the spin accumulation is directly 
calculated in an approximate but self-consistent manner. We also note that spin diffusion 



in a trilayer structure has been studied in a self-consistent manner in ref. [l7j, but only the 
linear instability was calculated (i.e. the threshold current for spin wave excitations), and 
the calculations were strictly for the case of a uniform current density flowing through an 
infinitely extended magnetic film without any consideration of boundary conditions. We are 
currently extending our work to the case of a trilayer, nanocontact structure so that the spin 
accumulation will be calculated in a self-consistent fashion, precluding the need to invoke a 
spin efficiency factor. 
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FIG. 7: (a) Nanocontact frequency spectrum as a function of dc applied current, (b-d) Local 
dispersion curves evaluated above, (solid blue curve), below (dashed black curve), and far from 
(dash-dotted red curve) the nanocontact along with numerically calculated wavenumbers and fre- 
quencies, (b) The circle is the wavenumber, frequency pair {k~,f~) calculated from the localized 
wave in Fig. [3^. (c) The diamond (square) corresponds to the wavenumber calculated above (below) 
the nanocontact k~ {k'^) in Fig. [Jja. The two frequencies /" and are peaks in the spectrum at 
/ = 28.5. (d) The diamond and square correspond to the wavenumber, frequency pairs as discussed 
for (c) but at the current / = 36 mA in Fig. Hh. 

In summary, we have used micromagnetics to predict the generation of a variety of re- 
sponses in single layer, nanocontact spin torque devices. A coUimated spin wave beam was 
observed over a range of currents and applied field angles, with the direction of the beam 
determined by the applied field azimuthal angle. The interplay of the applied field and 
the Oersted field act to form a spin wave corral, effectively trapping excitations under the 
nanocontact except along the direction where the oscillation frequency matches available 
propagation states. 



14 



The authors thank Bengt Fornberg and Keith Juhen for suggestions involving the nu- 
merical method used here. We also thank the University of Colorado at Boulder's Applied 
Mathematics Department for the generous use of its computational lab. 



* Electronic address: hoefer@boulder.nist.gov Contribution of the U.S. Government, not sub- 
ject to copyright. 

[1] S. L Kiselev, J. C. Sankey, I. N. Krivotorov, N. C. Emley, R. J. Schoelkopf, R. A. Buhrman, 
and D. C. Ralph, Nature 425, 380 (2003). W. H. Rippard, M. R. PufaU, S. Kaka, S. E. Russek, 
and T. J. Silva, Phys Rev Lett 92, 027201 (2004). 

[2] J. C. Slonczewski, J Magn Magn Mater 159, LI (1996). 

[3] M. L. PoUanski and P. W. Brouwer, Phys. Rev. Lett. 92, 026602 (2004). 

[4] M. D. Stiles, J. Xiao, and A. ZangwiU, Phys Rev B 69, 054408 (2004). 

[5] Y. Ji, C. L. Chien, and M. D. Stiles, Phys. Rev. Lett. 90, 106601 (2003). 

[6] B. Ozyilmaz, A. D. Kent, J. Z. Sun, M. J. Rooks, and R. H. Koch, Phys. Rev. Lett. 93, 176604 

(2004) . B. Ozyilmaz and A. D. Kent, Appl. Phys. Lett. 88, 162506 (2006). 

[7] S. Adam, M. L. Pohanski, and P. W. Brouwer, Phys. Rev. B 73, 024425 (2006). H. Xi, Y. 

Yang, J. Ouyang, Y. Shi, and K. Gao, Phys. Rev. B 75, 174411 (2007). 
[8] S. Kaka , M. R. Pufall, W. h. Rippard, T. J. Silva, S. E. Russek, and J. A. Katine, Nature 

437, 389 (2005). F. B. Mancoff, N. D. Rizzo, B. N. Engel, and S. Tehrani, Nature 437, 393 

(2005) . M. R. Pufall, W. H. Rippard, S. E. Russek, S. Kaka, and J. A. Katine, Phys. Rev. 
Lett. 97, 087206 (2006). 

[9] C. J. Garcia-Cervera and X. Wang, J. Comp. Phys. 224, 699 (2007). 
[10] R. E. Arias and D. L. Mills, Phys. Rev. B 75, 214404 (2007). 

[11] M. A. Hoefer, M. J. Ablowitz, B. Ilan, M. R. Pufah, and T. J. Silva, Phys. Rev. Lett. 95, 
267206 (2005). D. V. Berkov and N. L. Gorn, J. App. Phys. 99, 08Q701 (2006). G. Consolo, 
B. Azzerboni, G. Gerhart, G. A. Melkov, V. Tiberkevich, and A. N. Slavin, Phys Rev B 76, 
144410 (2007). 

[12] M. A. Hoefer, Ph.D. thesis. University of Colorado, Boulder (2006). 
[13] A. Slavin and V. Tiberkevich, Phys Rev Lett 95, 237201 (2005). 

[14] I. N. Krivorotov, D. V. Berkov, N. L. Gorn, N. C. Emley J. C. Sankey D. C. Ralph, and 

15 



R. A. Buhrman, Phys. Rev. B 76, 024418 (2007). 
[15] S. E. Russek, unpublished. 

[16] M. R. Pufall, W. H. Rippard, and T. J. Silva, Appl Phys Lett 83, 323 (2003). 

[17] A. Brataas, Y. Tserkovnyak, and G. E. W. Bauer, Phys. Rev. B 73, 014408 (2006). 

[18] B. Fornberg, SIAM J. Sci. Comput. 16, 1071 (1995). 

Appendix 

The inverse Fourier transform of eq. ([1]) gives 

^1= f [ Mr'A',r)[KL{R)+KniR)]r'dr'd<p', 

Jo Jo (A.l) 

R{r', (j)'] r, 0) = \/ r"^ ^ r'^ — 2rr' cos(</) — </>'), 

where lengths have been normahzed by the contact radius r^,, i? is the distance between 
the reference (r', 0') and source (r, 0) points, r is time, and we have assumed the thin film 
limit with |f = 0. The kernels Ki^ and are associated with the left and right interfaces, 
respectively. Their general expression is 

^^^^ Uk)kdk ^ ^ D/iwon) 

r Jq k coth(//t/r) + r/6' * ' 

K =[k^ + (rrf)2]V', d = n/lsf, I = I'/n, 

a =r*(±(522 + Wom\\)/{27rD). 

A plot of rK{r)/a is shown in Fig. [H]for the typical parameters considered in this work. 
The average magnetization direction over the contact is 

Ull 1 /"^'^ 

nil = nil = - / / u{r',(j)',z,t)r'dr'd(f)\ (A.2) 

Fill ttJo Jo 

which we use as the orientation of the longitudinal spin accumulation. The total spin 
accumulation is 

fh* = m*^ + [mz{0) + mz(t)]u\\, Uj_ = u — {u ■ u\\)u\\, (A. 3) 

where the longitudinal spin accumulation is the sum of the contributions from each interface 
z = and z = t. 
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FIG. 8: Lateral spin diffusion kernel rK{r)/a multiplied by r to remove the singular behavior 
and scaled by the constant a to compare the kernel contributions from each interface {rKL{r) / ai 
at z = 0, dashed curve and rKji{r)/aR at z = t, dash-dotted curve) and their combination 
{r[KL{r)/aL + Kii{r)/aR], solid curve). 

For a uniform applied field of the form = /i=K[sin(6'*) cos('?/'^,), sin(6'^,) sin('?/^^.), cos(6'*)], 
the dispersion relation for exchange spin waves in a thin film is 

(A.4) 

[qk^ + cos(6'e — 9*) — cos(26'e)], 

where 9e is the equilibrium magnetization polar angle satisfying /i* sin(^e ^^*) ~ | sin(26'e) = 
0. 

We briefly discuss the numerical method we have used to solve eq. ([2]). The polar coordi- 
nate system is a particularly efficient and accurate choice for nanocontact simulations. The 
discretization we use is non- uniform in radius ("inner" and "outer" grids) and uniform in 
angle 

ri = q{i - 1/2) = ]^{drin + drout){i - 1/2) + 

]-{-drin + drout)w\n[coii\i(^- ^-^tlZ^)/ cosh(-^)], 

^ XV xv 

cj)j = -7r + (j-l)d0, d4>=—, % = !,..., N, j = l,...,M, 

where i and w are parameters determining the location and width of the smooth change 
from the fine inner grid spacing dvi^ to the coarser outer grid spacing dvout- The advantage 
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of this discretization is that we can solve on a uniform computational grid {i — 1/2, (pj), but 
the physical grid {ri,(f)j) is clustered in and around the point contact where the interesting 
dynamics occur. The outer grid supports the propagation of spin waves of the appropriate 
wavelength away from the point contact. We choose a domain large enough < r < L, L ^ 
1 so that wave reflections off the boundary do not affect the solution inside the point contact. 
For comparison, our simulations on a circular domain with diameter 2Lr* = 2 ■ 100 ■ 40 nm 
= 8 /im involve MN = 32 ■ 502 = 16, 064 grid points. A nonuniform grid in the Cartesian 
coordinate system covering the same domain would involve AN"^ = 1 million gridpoints while 
a uniform grid would require 16 million gridpoints. 

The boundary condition for eq. ([2]) is the Neumann condition ^{r = L^cjy^r) = 0. We 
approximate radial derivatives using sixth order finite differences and the angular derivative 
u^^ is calculated using a pseudospectral, fast fourier transform method. An explicit Runge- 
Kutta 4"^ order time stepping method is used to advance equation ([2]) forward in time. 
To avoid severe time step restrictions due to the small grid spacing near the origin, we 
apply a smooth, radial grid dependent angular mask Qi(k) at every time step that filters out 



numerically induced small wavelengths near the origin ISj. The mask applied to positive 
angular mode k at the radial grid point rj takes the form 

9iik) = ^ + ^tanh ( ^\ , ^ ) , ki 



2 2 V / 2ri 

and is evenly extended to negative k values. 

Numerical parameters we use are: dvin = 0.05, dvout = 0.25, w = 10, {i = 126, M = 32, 
N = 342) for Figs. M and Ul {i = 273, M = 64, = 455) for Figs. [3] and H L = 60, ko = 3 
(wavenumber cutoff in angular mask at r = ri), and A A; = 1 (width of wavenumber cutoff 
in angular mask). We find no significant change in the results for more accurate grids and 
filtering parameters. 

The convolution m*j_{u±} in eq. (lA.ip is evaluated on the computational grid using Simp- 
son's rule in both the angular and radial directions. To deal with the removable singularity 
when i? = in eq. (lA.ip . we subtract off the small R behavior 



K{R) = af{R)/R ~ a{l/R + \n{R)/h + G) + o(l), i? < 1, 
where G is a constant. The resulting integral, at one interface, that we evaluate numerically 
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IS 




JQ 



R 



aAE{r^)u± + jU^ [ [ \n{R)r'dr', 
b Jo Jo 

u'j_ =u_L{r', (f)', r), u_L = u±{r, 0, r), 
where 4i?(r^) = J^^ r' / R dr' d(f)' and E is the complete elhptic integral of the second kind. 
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